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Abstract 

We employ relations between spacelike and timelike deep-inelastic processes in perturbative QCD 
to calculate the next-to-next-to-leading order (NNLO) contributions to the timelike quark-quark 
and gluon-gluon splitting functions for the evolution of flavour- singlet fragmentation distributions. 
We briefly address the end-point behaviour and the numerical size of these third-order corrections, 
and write down the second moments of all four timelike splitting functions. In the same manner 
we re-derive the NNLO result for the Higgs-boson decay rate into hadrons in the limit of a heavy 
top quark and five massless flavours, and confirm the recent N^LO computation of this quantity. 



In this letter we present new results on the scale dependence (evolution) of the parton fragmentation 



distributions D^Ux.Q^). Here x denotes the fraction of the momentum of the final-state parton 



/ carried by the outgoing hadron h, and is a timelike hard scale, such as the squared four- 
momentum of the gauge boson 'm ^ y,Z ^ h+X. This scale dependence is given by 



where the summation over j — q, q, g is understood. The timelike splitting functions Pjf admit 
an expansion in powers of the strong coupling ttg, 

Pf, (x,a,(e2)) = a,pf\x) + alpf\x) -f alpf\x) -f ... , (2) 

where we normalize the expansion parameter as = CXs(2^) / (47r). The leading-order (LO) terms 
in Eq. Q are identical to the spacelike case of the initial-state parton distributions, a fact often 

referred to as the Gribov-Lipatov relation [1]. Also the next- to-leading order (NLO) contributions 

(1) r 

Pji {x) have been known for more than 25 years. These quantities are related to their spacelike 
counterparts by a suitable analytic continuation [2-4], see also Ref. [5]. 

In a previous publication [6] we have calculated the next-to-next-to-leading order (NNLO) 
splitting functions Pns {^) for the non-singlet combinations of quark fragmentation distributions. 
This calculation was based on an analytic continuation of the corresponding unrenormalized par- 
tonic structure function in deep-inelastic scattering (DIS) [7], performed after subtracting contri- 
butions due to the quark form factor [8] and backed up by another relation between the spacelike 
and timelike cases conjectured in Ref. [9]. We now address the flavour-singlet timelike evolution 



Here ® abbreviates the Mellin convolution written out in Eq. ([T]), and rij: stands for the number 
of effectively massless quark flavours. Specifically, we extend the approach of Ref. [6] to derive 
the NNLO diagonal entries P^^^ {x) and Pg^ {x) in Eq. ([3]). The relation of this calculation to 
Higgs-boson decay will be addressed below. 

For brevity focusing on the gluonic case, our calculation starts form the unrenormalized struc- 
ture function F^^ for DIS by the exchange of a scalar (|) coupling (like the Higgs-boson in the 

heavy-top limit) directly only to gluons via ^C^^Gc^, where G^y denotes the gluon field strength 
tensor. This quantity has been computed to three loops for the determination of the spacelike 
NNLO quark-gluon and gluon-gluon splitting functions in Ref. [10]. In dimensional regulariza- 
tion with D = 4 — 2£ its perturbative expansion in terms of the bare reduced coupling can be 
written as 



F^_^{alQ^) = 5(1 -x) + £ {a^r ^ " 

1=1 VA* / 



The n-th order terms F^g^^ are then iteratively decomposed into contributions arising from the 
analogous expansion coefficients !F,„<n of the ^gg form factor [1 1] and remaining 'real' parts 2?,„, 
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F^^l'^ = 2iFi5(l-x)+^i 

= 2j25(l-x) + (^i)26(l-x) + 2iFi^i+^2 

= 2535(l-x)+2Ji525(l-x) + (2iF2 + (Ji)')^i+2Ji^2 + ^3 • (5) 

Note that the functions 3(,,„>2 do not only collect tree-level amplitudes but also combinations of 
real-emission and virtual corrections. As in Ref. [6], this will lead to a problem in the third- 

(2)T 

order analytic continuation. For igg , however, this problem can be fixed afterwards in complete 
analogy to the previous non-singlet quark case (see below), thus we can ignore it for the moment. 

The analytic continuation of the form factor to the time-like case is known. The x-dependent 
functions ^„ are continued from x to l/x [2-5], taking into account the (complex) continuation of 
(see Eq. (4. 1 ) of Ref. [8] ) and the additional prefactor x^^^^ originating from the phase space of 
the detected parton in the timelike case [12]. We have performed this continuation using routines 
for the harmonic polylogarithms (HPLs) [13] implemented in FORM [14]. The only subtle point 
in the analytic continuations is the treatment of logarithmic singularities for x ^ 1 starting with 

ln(l— jc) \n{l — x) —lnx + in . 

After these analytic continuations the one-parton inclusive fragmentation function in (|)-decay is 
re-assembled order by order analogous to Eq. ([5]), keeping the real parts of the continued ^„ only. 
Then the renormalization of the operator G^yG^^ and the strong coupling constant is performed. 
Finally the timelike splitting functions (and coefficient functions) can extracted iteratively from the 
mass factorization relations 

^<s?,8 - "e^gg +%g +^ %g +■•• 

+ e|.<f-/.f*<fU... (7) 



6£3 



6£2 

_ 1 /2p(2)^ , 3 pWT WT . p{0) {2)T _ (0) / (0) 3^5..\ 

(3)r 1 (1)7- (i)r p(0) (2)r 1 (0)/ (0) „ . . 

Here all products of x-dependent (generalized) functions are to be read as Mellin convolutions or 
as products in Mellin-A'^ space, employing routines for harmonic sums and their inverse Mellin 
transform back to .x:-space [14, 15]. Obviously the determination of igg from Eq. (|8]) requires 

{2)T (2)T 

the 'off-diagonal' two-loop coefficient function ^ . This quantity, ^ and the corresponding 
first-order functions (for all these our normalization differs from the standard convention by a 
factor of two) can be calculated via a direct analytic continuation of ^,|,''^^'^'*- We have checked this 
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fact by comparing the corresponding gluonic results for the photon-exchange case to an explicit 
two-loop calculation to order £ [16]. At the third order in tts, however, this direct continuation 
fails to correctly reproduce the contributions already at order e^, thus we cannot derive the 
off-diagonal timelike NNLO splitting functions in this simple manner. 

We are now ready to present our results for the diagonal timelike splitting functions. For 
completeness we start at NLO where we, of course, reproduce the results of Ref. [3]. Adopting the 
notations of Ref. [13] for the HPLs, the timelike - spacelike differences can be written as 



- 20/9 jc"^ -3-X + 56/9 - {3 + 7x + S/3x^) Hq + 2(1 +x) Ho,o 



(9) 



8Ci(^;7gg(x) [ll/3 Ho -4(Ho,o + Hi,o + H2)J +[6(1 -x) -22/3(^-1 -x2)]Ho 
-8(l+x)Ho,o) - 16/3C4n^pgg(jc)Ho + SCpnf(20/9x-^ +3+x-56/9x^ 
+ [4 + 6x + 4/3 (jc- 1 + ] Ho + 2 ( 1 ;c) Ho,o^ 



(10) 



where we have used the abbreviation 

PggW = 'i-fi'i- —x) + 1/x — 2 + X — x^ . 
Note that the non-HPL terms in Eqs. Q and (flOl) are identical up to an overall sign, cf. Ref. [9]. 

T S 

The functions Pps denote the 'pure singlet' contributions from which the quark-quark entries in, 
e.g., Eq. (|3]) are obtained by adding the corresponding non-singlet quantities. 

The difference between the timelike NNLO pure-singlet splitting function and its spacelike 
counterpart of Ref. [10] reads 



+ SCj^Cpnf (^269/6jc"1 + 14 + 1 13/2x - 346/3 + ^2 (172 + \67x + Sx^)/3 
-^3(12x"i-13 + 65x-28x2)-2(l+jc)[l6^|-F4H_i,o,o + 9H3,o + 4H3j 
+ 1 C2H2 - 1 2 H2,o,o - 2 H2, 1 ,0 - 6 H2,2 - H4] + 8/3 (x- 1 + x^) [4 H_ 1 ,0,0 + C2H0 
-2(1 -x) [8(H_3,o + H_2Ao)+5^2Hi-9C2Ho + 25/12Hi,o-6Hi,o,o-Hli,o 
- 3Hi,2 + H2,i] + S/3{x-^ -x^) [6Hi,o,o + Hi, 1,0 + 3Hi,2 - 5 ^2Hi - H2,o - H2,i 
+ 2/3 (4x^1 + 27 - 63x + 28^2) H_2,o - 2/3 (20jc"^ - 21 + 9x + 56x^) H_i,o 
+ (89/9x"^+55 + 1021/6x + 2297/18jc2-46^3-22x^3)Ho-(8/9jc"^ 
+ 293/6 + 370/3^ + 538/9x2 + 2^2(1 - 7x)) Ho,o - 32xHo,oao + (5 - \6/3x-^ 
+ 85x) Ho,o,o - 1/6 (1 15x^1 + 362 - 292x - 185 jc^) Hi - (6x"^ +48 + 59x 
+ 22x^) H2 - 2 (5 + X - 8/3 jc^) H2,o + 4 (2/3 jc" ^ + jc + 2^2) H3) 



+ 8C^n^(^ -217/18- 55/3x^1 + 122/9X+ 101/6x2 + ^3 (16x"^ +36 + 24x) 

- ^2 ( 127 + 1 88 X + 1 28 jc^) /3 + 2 ( 1 + x) [l 6 ^1 + 1 7 ^3Ho + 8 ^2 -^Ho - 7 ^2Ho,o 
+ 10 C2H2 + 9 H2,o - 12 H2,o,o - 2 H2, 1,0 - 6 H2,2 + 9 H3,o + 4 H3, 1 - H4 
+ 2(1 -x) [5 C2H1 + 2H2,o + 139/12Hi,o - 6Hi,o,o - Hi,i,o - 3Hi,2 + H2,i 
+ 8/3 {x-^ - x^) [5 ^2Hi + 5/3Hi,o - 6Hi,o,o - Hi,i,o - 3 Hi,2 + 2H2,o + ^2,1 

- (527 + 2473 x + 8 1 1 jc V 72 ^2) / 1 8 Ho + (62 + 8 1 /2 jc + 208/9 jc^) Hq^o 
+ (6 + 18x- 8^2) Ho,o,o + (385/18x^1 + 190/3 - 143/3X- 667/18^2) Hi 
+ (28/9x"^ +71 + 46x + 248/9 jc^) H2 - 4/3(4x"' - 6 - 3x + 8jc2)H3) 

+ 8Cpn} (2/9 (23x - 2x-i - 20 + 2(1 +x) [^3 - ^2Ho - Hi +H2,o + H3 

- Ho,o,o] - (1 - x) (Hi - Hi,o) + 4/3 {x'^ -x^) Hi,o + 2/9 (3 + 18x + lOx^) Ho 
-(7+x-4x2)/3Ho,o-(20x-i-56x2)/9Hi + (3 + 7x + 8/3x2)(^2-H2) 

The corresponding result for the gluon-gluon splitting functions is given by 



hP^^ix) 



+ 16Ci (^pgg(x) [(1025/54 - 1 1/3 C2 - 2^3) Ho - 49/3 Ho,o - 33Ho,o,o + 16Ho,o,o,o 

- (268/9-8^2) (Hlo+ H2) -44/3 (Hi,o,o + H2,o + H3) + 12Hi,o,o,o + 4H2,o,o 
+ 4H3,o + I2H4] + Pgg{-x) [l6H_3,o - 16 C2H-2 - 16H_2,-i,o - 22/3 H_2,o 
+ 28H_2,o,o + 8H_2,2-16H_i,_2,o-32H_i,_i,o,o-16C2H-i,o-44/3H_i,o,o 
+ 36H_i,o,o,o + 8H_i,2,o + 16H_i,3 + (14^3-ll/3C2)Ho+16C2Ho,o + llHo,o,o 

- 16Ho,o,o,o - 4H3,o - I2H4] + (1 +x) [ - 24H_2,o - 48H_i,o,o + 14/3 H2,o 

+ 28/3 H3I + ( 1 - x) [32 (H_3,o + H-2,0,0) -(881/36- 24^3) Ho - 27 (Hi ,0 + H2) 



-44/3 {x-^+x^) 



+ ix-'-x^] 



226 1/54 Ho 



2H_2,o + 4H_i^o,o — H2,o — 2H3 
+ 134/9 (Hi.o + H2)j - (44x^1 + 86 + 14x + 132jc2)/3 ^2Ho + (536jc"^ + 425 
+ 515x + 752x2 + 288C2)/9Ho,o + (88x"^-10 + 8x + 44jc2)Ho,o,o + 64xHo,o,o,o 
+ 16Q\(pgg(x)[- (158/27 -2/3^2)Ho -4/9 Ho,o + 6Ho,o,o + 40/9(Hi,o + H2) 
+ 8/3 (Hi^o,o + H2,o + H3)1 + 2/3 Pgg(-x) [2H_2,o + 4H_i,o,o + C2H0 - 3Ho,o,o 



(1+x) ^2Ho-H2,o-2H3 -(1-x) 173/9Ho + 2(Hi^o + H2) +{x 



-1 -x^] 



913/54 Ho + 26/9 (Hi,o + H2) + 4/9 (35 jc^^ + 21 + 48 jc) Ho,o + 4 (1 + 4jc) Ho,o,o 



16 



+ —CAnf[p^^{x) 



10Ho + 12Ho,o 



+ 12(l+x)Ho,o + (13(x-i-x^)-9 + 9x)Ho 



(11) 



+ 8QC^n^ - 2/?gg(x) Ho - 269/6x"i - 14 - 1 13/2x + 346/3^2 - ^2 (172 + 167x 
+ 8x2)/3 + ^3(12x"^-13 + 65jc-28jc2)+2(1+x)[i6^| + 2 H_2,o - 4 H_ 1 ,0,0 
+ 17 C3H0 + 4/3 C2H0 - 3 ^jHo.o + IOC2H2 - 12H2,o,o - 2H2,i,o - 6H2,2 + 4H3,i 
+ 9H3,o-H4] +8/3(x-i+x2)[4H_i,o,o + C2Ho] -2(1 -x) [8(H_3,o + H_2,o,o) 

+ 18 H_2,o + 9 C2H0 + 6 C3H0 + 4 C2Ho,o - 145/ 12 Hi,o] + [ ^ (x"' - x^) + 2( 1 - x) 
3H_2,o - 1 1/3 Hi,o + 5 ^2Hi - 6Hi,o,o - Hi, 1,0 - 3Hi,2 +H2,i] + (40^-1 - 54 
+ 18x+ 1 12x2)/3H_i,o - (59x^1+45 + 1081/6x+ 157/2^2) Hq - (464x^1 
+ 329/2 - 146x - 66jc-)/9Ho,o - (80/3x"^ - 17 + 15x) Ho,o,o - 32xHo,oao 
+ (1 15jc"1 + 362 - 292jc - 185x2)/6Hi - 1/9 (34jc"1 - 546 - 417jc - 286^2) H2 
+ (8x^1 + 10 - 14x - 24x2)/3 H2,o - 8/3 (x^^ + 5 + 13/2x + 3x2) j 

+ 8C|n^ (^217/18 + 55/3x^1 - 122/9x - 101/6x2 - ^3 (16x"i +36 + 24x) +^2/3 
(127 + 188x + 128x2) - 2(1 +x) [16 + ^2Ho,o + 10^2H2 + 17 ^3Ho - 12H2,o,o 

J(x-i-x2)+2(l-x) 5^2Hi 



+ (4x"^ + 283/6 + 239/2x + 739/18x" 



- 2H2,i,o - 6H2,2 + 9H3,o + 4H3,i - H4 

+ 3Hi,o-6Hi,o,o-3Hi,2-Hij,o + H2,i 

- 8 ^2 - 20^2-^ - 16/3 ^2^^) Ho - (18 + 97/2x+ 16x2) Hq^q - (6 - 6x + 8x2) ^^^^ ^ 

- (385 x"^ + 1 140 - 858x - 667x2)/18 Hi + 53/6 (1 -x)Hi,o - (20/3x"i +45 
+ 72x + 24x2) H2- (32/3x^1 + 14 + 6x) H2,o - (16/3x"i - 8 - 12x) H3) 

+ 8/9 Cpnf (4x-^ + 40 - 46x + 2x2 - 9^2 (3 + 7x + 8/3 x2) - 6( 1 + x) 3 ^3 + ^^^^ 

(92/3x"^ - 6 + 48x- 32/3x2) Ho -(16x^^+83 



+ 3Hono-H2n-5H3 



+ 101 X + 28x^) Ho,o + (20x^1 + 27 + 9x - 56x^) Hi + (4x"^ + 3 - 3x - 4x^) Hi,o 



-1 



+ (16x"^+39 + 51x + 8x^)H2) . 



(12) 



Eqs. (fTTI) and (fT2l) represent the main new results of this letter. Analogous to the non-singlet result 
of Ref. [6], the coefficient of C| pgg(x)//o,oC2 in Eq. (fT2l) differs from the result of the analytic 
continuation as specified above. We have determined the correct coefficient via the momentum 
sum rule for njr = 0. The complete C|, C^n^- and C^nj parts of Eq. (fT2l) have been independently 

derived by applying the approach of Ref. [9] to the (non-singlet like) non-C^ contributions to igg. 

(2) (2) 

Furthermore an inspection of the above results reveals that also 5igg + 5ips does not receive any 
non-HPL contributions. Since the determination of these two functions uses largely independent 
information, we view this fact as another non-trivial check of our calculations. 

(2)T 

As predicted in Ref. [9], the large-x behaviour of igg is identical to that of its spacelike 
counterpart (see Eq. (4.16) of Ref. [10]) up to the sign of the subleading ln(l — x) contribution. 
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The pure-singlet splitting functions are suppressed by two powers of {I —x) for jc ^ 1. We now 
turn to the small-x limits. Up to terms suppressed by powers of x the NLO quantities read 

80 



XPc 



qq ^ ('"•) 



16CjLo - 



8 



1 1 Q + 2ny (Q - 2C^)J Lq - — (Q - 2Cf 



92 



(13) 



with Lq = \nx. The log-enhanced contributions of our new diagonal NNLO splitting functions are 

(14) 



xPqq {x) 



32 „„ .--3,. 2. 8 



xP^'^'(x) 



^ Cj,Cpnf{2Li + Lt) + -(155 + 72^2)QC^n^Lo + 0(i; 
^ Ci l4 + ^ (33 Ci + 6Ci/7^ - lOQC^n^) 4 



3 
8 

+ 9L 
8 

+ 27 



(389- 144^2) Cj + Ueclrif - ISlC^CpHf + AnjiC^-lCp 



(4076 -990 ^2 -972 ^3) + (739-36^2)Ci7i^ 

2„ I /tz:„2 



- (1819- 144^2) QC^n^ + 108C;/i^ + 46/i;(Q - 2C^ 



\f 



Lq + 0(1) . (15) 



Thus, in contrast to the corresponding spacelike quantities, the NNLO timelike singlet splitting 
functions receive double-logarithmic contributions with very large coefficients, with the leading 
term of Eq. (fT5l) agreeing with Refs. [17, 18]. Consequently both xP^^ and xPg^ - despite a 
large cancellation between the leading- and subleading logarithms for the latter quantity - show 
a huge enhancement already at x > 10^ as illustrated in Fig. [T] For this and the next figure the 
harmonic poly logarithms have been evaluated using the FORTRAN package of Ref. [19]. 

Returning to the region of medium and large values of x. Fig. |2] shows the LO, NLO and 
NNLO approximations to the diagonal entries in Eq. ([3]) at a scale relevant to gauge-boson and, 
maybe, Higgs decay. Also here the higher-order corrections are larger than for the spacelike case, 
in particular for the gluon-gluon splitting. Nevertheless the perturbative expansion appears well- 
behaved at least for .x: > 0.1. As usual, the region of safe applicability of the NNLO approximation 
for Eq. ([D) will be wider as an effect of the Mellin convolution. 

As mentioned above, we are presently not in a position to derive the ^2-terms of the off- 
diagonal quantities iqg {x) and /gq {x) (with the exception of some /ly^-enhanced contributions). 
However, their non-^2 second moments provide another check, via the momentum sum rule, of 
our new results ([TT]) and ([T2l) . Furthermore, using these results and the momentum sum rule, the 
missing off-diagonal terms can be reconstructed at = 2. Hence we can write down the complete 
NNLO expressions for the second moments of all four splitting splitting functions. 



qq 



{N = 2) 



gq 



{N = 2) 



r3 



/ 54556 7264 



(- 



2^ /^6608 



\ 243 



55 



2432 ^ 2464 ^ 
-9- ^2 + -9- ^3 

296 . 512. 



243 
128 



27 



^2-320^3 + 256^1 



Cpflj: 



.(10910 6A 



2281 32 64 
^-y^2 + yC3 



(16) 
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In'^x 



* 1/2000 



_l I I I I I I I I 



-20 

10'^ 10'^ 1 10'^ 10'^ 
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X X 
Figure 1: The third-order timelike quark-quark and gluon-gluon splitting functions for five 
flavours, multiplied by x ( 1 —x) and divided by 2000 ~ (47r)^ for display purposes. Also shown are 
the respective leading small-x contributions and the corresponding spacelike splitting functions. 
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Figure 2: The perturbative expansion of the timelike quark-quark and gluon-gluon splitting func- 
tions, again multiplied by jc(l — jc), at a typical value of the strong coupling constant. 
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^2/^2 160^ 64 ^\ /2681 760 ^ 56. 



. ,10570 352 ^ 32^, ^ 2/^41 128 ^ \ 

After inserting the numerical values of the QCD colour factors and the Riemann ^-function, these 
results yield, for five light flavours, the benign perturbative expansion 

Pg^(A^=2, nf = 5) ~ 8as/(97r) (l - 0.687 a,, + 0.447 + ...) 

Pq^(A^=2, n^ = 5) ~ 5as/(67r) (1 - 1.049as + 1.163a2 + ...) . (18) 

We finally turn to the timelike gluon coefficient functions in Eqs. © - ([8]) and their quark 
analogues c|'q^(A:). In the limit of a heavy top quark and negligible masses of all other flavours, 
the coefficient functions for the standard-model Higgs-boson differ from those of this scalar (|) only 
by a perturbative prefactor known to N-^LO [20]. Thus the latter quantities are directly relevant to 
one-particle inclusive Higgs decay, and the second-moment combination 

{Clq + Cl^){N = 2) = \+aJ^^ +alcf^ +alcf^ + ... , 

in] 

directly enters the Higgs decay rate in the above limit. For the expansion coefficients our 
analytic continuations lead to 

fi 



73 „ 14 
T 



= -CA--nf. (19) 



(2) ^2/^37631 242 . ^ /6665 88 ^ 

= — C2-110C3j-Q.,^^--C2+4C3 



131 \ , ^5_08 _ 8 

27"~ 3 



C^nW— -24^3 --^2) , (20) 



(3) , 15420961 178156 ^ 3080 ^ \ ^ 2/^413308 56 ^ 
= M2)+Cl^-^^ 2^^3 + ^C5j + Qn2(^^^ + -C3 

. /2670508 9772 ^ 80 ^ \ /23221 ^^^^^ ^^^^ 
- i—^ ^ ^3 + Y ^5 J - CpC^Hf 1364 ^3 - 160 ^5 

+ - 320 + 192 C3) + C^n} (440 - 240 ^3) - n'f - g ^3) (21) 

with a function /(^2) which we cannot derive at this point. Eq. (fT9l) reproduces a well-known NLO 
result of Refs. [21,22]. Our coefficient (l20l) represent a completely independent re-derivation of 
the NNLO expression first obtained in Ref. [23] (see also Ref. [24]). Finally the third-order result 
([2T)) provides, despite the missing ^2 -contributions, a highly non-trivial confirmation of the recent 
N^LO calculation of Ref. [25] to which the reader is referred for a further discussion. 
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To summarize, we have used relations between spacelike and timelike quantities to compute 
some important higher-order QCD corrections for timelike processes. We expect that further 
progress can be achieved along these lines. However this will require improvements on the present 
decompositions into purely real and real-virtual terms which are beyond the scope of this letter. 

Form and Fortran files of our results can be obtained from |http: //ar Xiv.orgl by down- 
loading the source of this article. Furthermore they are available from the authors upon request. 
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